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ABSTRACT 



Aims. We investigate the chance of detecting proto-planetary or debris disks in stars that induce microlensing event (lenses). The 
modification of the light curves shapes due to occultation and extinction by the disks as well as the additional gravitational deflection 
caused by the additional mass is considered. 

Methods. The magnification of gravitational microlensing events is calculated using the ray shooting method. The occultation is taken 
into account by neglecting or weighting the images on the lens plane according to a transmission map of the corresponding disk for a 
point source point lens (PSPL) model. The estimated frequency of events is obtained by taking the possible inclinations and optical 
depths of the disk into account. 

Results. We conclude that gravitational microlensing can be used, in principle, as a tool for detecting debris disks beyond 1 kpc, but 
estimate that each year of the order of 1 debris disk is expected for lens stars of F, G, or K spectral type and of the order of 10 debris 
disks might have shown signatures in existing datasets. 

Key words. (Stars:) circumstellar matter - Gravitational microlensing - Astrometry 



c3 1. Introduction 



Since the discovery of the first circumstellar dust disk around 
Vega by the Infrared Astronomy Satellite (Aumann et al j 1984b . 
at least 101 circumstellar disks have been resolved around pre- 
main and main sequence stars out to a distance of 1 kpcQ in- 
cluding 17 so-called debris disks. Proto-planetary disks may be 
quite common during the pre-main sequence life of a star, but 
with a maximum life span of only 10 Myr, a galactic star forma- 
tion rate of 5 stars /yr, and an upper limit of 400 billion stars in 
the MilkyWay (McKee & Williams [19971 Diehl et al. 12006b . the 
chance of a random star containing a proto-planetary disk is only 
about 0.01 %. The debris disks that remain thereafter last much 
longer, perhaps as long as 10 Gyr (e.g. Greaves 2005 ), and so 
should be much more common. Trilling et al. (2008 ) have shown 
that the occurance of 70 fim excess emission in a sample of F, G, 
and K stars - presumedly due to dust in some form of debris disk 
- is 16 ± 3%. In contrast, practically no main sequence M-stars 
show excess 70 [im emission: Rhee et al. d20071 > have found only 
one disk (AU Mic) from a sample of ~ 900 M-stars. 

While a common phenomenon, circumstellar disks around 
main sequence stars can be difficult to detect using standard 
methods, coronography and the detection of excess infrared 
emission. A third method - detecting the disks in absorption 
against a background source - works very well for dense proto- 
planetary disks seen against an emission nebula ("proplyds"; 
McCaughrean & O'Dell 1996) but debris disks around main se- 
quence stars are optically thin and unlikely to be found in front 
of a bright background. 

Gravitational lensing could provide an additional method for 
detecting and characterizing circumstellar disks by creating a 
well-defined and bright if temporary and geometrically compli- 
cated background source. If the projected distance between the 
lens star and the background source is small enough, the source 



star appears to separate into two geometrically extended objects 
- the source of the photometric magnification in so-called mi- 
crolensing, where the lensed object is not resolved. The angular 
scale of gravitational lensing is given by the Einstein radius 



4GM D LS 
c 2 D S D L 



(1) 
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where Dls denotes the lens-source distance, Di the lens- 
observer distance, D s the source-observer distance and M is the 
lens mass. For a typical galactic microlensing event observed to- 
wards the galactic bulge (Paczyriski 1991) with a 0.5 M lens 
at 6 kpc and a source at 8 kpc, one obtains an angular scale of 
the order of a mas (Wambsganss 2006), corresponding to a dis- 
tance of ~ 2.5 AU at the distance of the lens. Thus, the typical 
microlensing events studied primarily to detect low-mass exo- 
planets around lens stars (Beaulieu et al. 2006; Bennett et al. 
2008) may also contain information about much more diffuse 
circumstellar matter as well. 

Bozza & Mancini (2002a) and Bozza et al. (|2002b) have 
studied the effects of microlensing by lenses consisting of gas 
clouds and interacting binaries which may be used as an approx- 
imation for a face-on proto-planetary disk and Zheng & Menard 
( 2005) have considered the effects of disks around microlensing 
source stars, but to our knowledge the detectibility of circum- 
stellar disks around lens stars has not been considered in the lit- 
erature. About 70% of the observed lens stars towards Baade's 
window are M stars and smaller and only ~ 25% are of spec- 
tral type F, G, or K (Dominik 2006). Given that the microlensing 
survey groups OGLE (Udalski et al. 1 19941 ) and MOA (Abe et al. 
119971 ) are detecting about 1000 events per year, the fraction of 
events containing a debris disk should be significant: about 4.0% 
or currently ~ 40/year. The chance of detecting proto-planetary 
disks among lens stars is correspondingly lower: 0.0025% or 
~ 0.25 per year (neglecting the additional absorption likely to 
be found within star-forming regions). Although only 1 in 8000 
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Fig. 1. Light curve of the point source point lens model occulted by an optically thick edge-on disk is plotted along with the 
corresponding source and image tracks and the astrometrically detectable astrometric offset from the source. Features A and B are 
caused by the occultation of the less magnified image while feature C is obtained from the occultation of both images. The impact 
parameter is 0.4 #g, the inclination 89.8 ° and the disk is rotated by 15 °. 



lensed stars have proto-planetary disks, future surveys like the 
Microlensing Planet Finder (Bennett et al. 120041 ) or EUCLID 
(e.g. Beaulieu et al.2008'1 have a chance of detecting a signifi- 
cant number of such rare events. 

The effects of a circumstellar disk on microlensing 
lightcurves can be separated into two very different regimes. 
The disk can occult or at least extinct the light of the source 
star as it is bent around the lens star - a simple geometric effect 
given a description of the geometry and matter distribution of 
the disk and the geodesic paths of the photons from the source. 
Alternatively, the mass of the disk itself can distort the geodesic 
paths of the photons and hence change the apparent areas of the 
distorted source images and the resulting lightcurve independent 
of any occultation or extinction effects. The purely geometric 
disk parameters are the projected size, shape, inclination, and 
orientation angle <f> of the major axis of the disk relative to the 
path of the source star and the size of the lens star's Einstein ra- 
dius. The typical outer radii of proto-planetary disks are about 
100 - 8 00 AU and of debris disks about 50 to 100AU, (e.g. 
Greaves 2005) corresponding to 20 - 320 Be for a typical lens. 
In both cases we adapt an outer radius of 150 AU as chosen by 
Brauer et al. (|2008t . 

In the following we will discuss the chance of detecting cir- 
cumstellar disks by simulating a variety of potential light curves. 
We will start with optically thick disks, where the optical depth is 
greater than one (cf. Rrivov 2005) representing the young proto- 
planetary disk phase and show that the corresponding magnifi- 
cation structure is scaled according to the ratio of the disk mass 
and its stellar host. This insight leads us to the most likely case of 
a microlens surrounded by an old debris disk modelled as opti- 
cally thin and geometrically thick attenuating disks with optical 
depths t ~ 0.01. 



2. Optically thick proto-planetary disks 

Given that the major microlensing effects occur for impact pa- 
rameters less than a few Einstein radii, a magnified background 
source will only be visible behind proto-planetary disks with 
high inclinations; disks with low inclinations may be too large 
to show any relativistic effects and will only show simpler oc- 
cultations of non-magnified sources. Similarly, debris disks with 
high inclinations are more likely to be detected because both the 



projected extinction and the projected mass-density responsible 
for any additional lensing effect will be higher. 

An optically thick disk can be modeled as an occulting disk 
with given inclination and orientation angles as parameters. The 
features of the corresponding light curves and the astrometri- 
cally observable change of the centroid of both images are il- 
lustrated in Fig.Q] neglecting the gravitational deflection due to 
the disk mass. The less magnified image is occulted close to the 
maximum of magnification. Because the disk structure is bound 
to the lens, it is unlikely to occult the magnified image exclu- 
sively; unless the edge-on-disk is nearly parallel to the straight 
line through the two images and the lens. An asymmetric total 
occultation of both images as in Fig.[T]indicates that an optically 
thick disk structure is present. 

In addition to the purely photometric effect, the occultation 
also produces a shift in the center-of-light position of the lensed 
source, albeit small. The astrometric offset due to the occulta- 
tion of the less magnified images, without parallax effects (see 
Gould 119921 ), can be estimated by assuming that these images 
are located close to the Einstein radius 



«_ 1 ( u Vm 2 + 4^1 

60 * —Op = - 1 

H h 2{ u 2 +2 j 



where u denotes the lens-source separation in /i the total 
magnification and fi ± the magnification factor of magnified and 
demagnified image. For a typical galactic microlensing event, 
the angular Einstein radius is w 0.4 mas: if the source-lens sep- 
aration is below 1 Einstein radius, the pre-factor in Eq. [2] takes 
values between 0.12 and 0.5 and the total offset is of the order 
of 10/ias Fig. [TJ, which is detectable using the planned GAIA 
(Lindgren & Perry man 1996) or SIM (Unwin et al. E(X)8]) satel- 
lites. For ground based observations Rattenbury & Mao (2006) 
propose the usage of closure phase interferometry which could 
provide the desired accuracy in this context. 

Given the rareness of such young disks and the difficulty 
of seeing microlensing effects, detecting such disks photomet- 
rically and astrometrically is not impossible but unlikely. 
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3. Optically thin debris disks 

3. 1 . Geometrically thin disks 

Up to now, we have neglected the gravitational deflection by the 
mass content of the circumstellar disk, represented by the mass 
surface density 



s \ 0.001/ 0.5mJ\ 



Rr 



-2 



150 AU I 



r 

R~7) 



-0.8 



(3) 



(e.g. Brauer et al. 12008 ), where qo = Mo/M+ is the ratio of the 
disk mass to the stellar mass, and Ro is the effective outer radius. 

Assuming qo is small, the maximum change in deflection 
angle and magnification is for edge-on oriented disks, whose ef- 
fective linear mass distribution can be calculated by integrating 
the surface density. The deflection angle of the combined point 
source point lens (PSPL) and edge-on disk model has no rota- 
tional symmetry, so we have to consider the vectorial deflection 
angle a. An analytical expression can be derived if one integrates 
the normalized line density in the lens plane A(9i ) along the disk 
radius 9o, expressed in units of the Einstein radius (cf. Narayan 
& Bartelmann [T996t : 



a{0) 



1 



(l + qa> 



|0| qo + 1 J-, 



wo 

-ft, |0-0 



1/2 l 



(4) 



The normalized line density for a mass distribution parallel to 9\ 
is given by 



(5) 



where 9\z define the lens plane and the integral over #3 projects 
the surface density distribution onto the line. The resulting ex- 
pression contains the hypergeometric function 2F1 (Abramowitz 
& Stegun [T972l : 



W) = I0i 1 1/5 



jr 1 



0.45(22.18(6^ -9 2 f /l °9 



1 



1 4 9 



9\ ^ 



18.94 



(6) 



10' 10' 10' 9 2 D -9\ ) 
which can be efficiently evaluated sufficiently away from the 
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Fig. 2. The normalized line density of an edge-on oriented disk 



with a surface density E oc y-^jj) > ™ a disk with r, 
0.03 AU and 9 D = 60 % is plotted. 



in 



singularities at ±9d- When a third order polynomial fit is used 
outside a 40 0g radius, as shown in Figf2] the total deflection 
angle (Eq. |4|i can be numerically integrated. 

Even though the lens equation with the additional term for 
the deflection angle cannot be solved analytically, the disk can 
be embedded into a ray shooting framework (Kayser et al. 1986 ). 
For qo ~ 0.2, the calculated magnification offset within 1 9g 
from the disk is ~ 1 and hence potentially observable. Close 
to the linear mass distribution representing the disk, one sees 
an additional rectangular magnification caused by the attraction 
of rays (Fig. |3). For small source-lens separations u (here u < 
0.05 9e) the singular peak splits in two at the central region. For 
qo ~ 0.20, the affected region of the extended box is 0.02 9e and 
the central distortion has a relative deviation larger than 1 % for 
0.2 9e, where the Einstein radius is calculated for the total mass 
of the system. 
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Fig. 3. Top: magnification map for a linear mass distribution cre- 
ated by ray shooting (Eq. |4]) for the combination of a point lens 
and a disk mass of 20 % of the stellar host along with contour 
lines of magnifications at fi = 2.5, 3.7, 4.8 and 6. Profiles ex- 
tracted from the magnification map located in the source plane 
correspond to light curves. Middle: change relative to a point 
source point lens (PSPL) model. Bottom: vertical cuts corre- 
sponding to deviations in the light curves for two perpendicular 
trajectories of the lens and source for 0, 0.5. &g 



According to Eq.[4] the width of the distortion is proportional 
to A(9i) and thus to the mass ratio qo for ijo « 1. The rectan- 
gular distortion is redistributed in case of an inclined disk due 
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to the conservation of rays: for i - 89° the conserved integral of 
the rectangular distortion reduces the peak deviation to 2% and a 
face-on rotational symmetric disk has no effect at all. If the disk 
has no central gap, a deviation containing a splitted singularity 
is observable as long as the finite size of the source star does not 
smooth it out. Considering a bright giant source star with 10 R 
and thus 0.014 6e at 8 kpc, this limit is reached at a mass ratio of 
0.014 and, in case of inclined disks, even earlier. 

The calculation of the extinction in a geometrically and op- 
tically thin debris disk requires an estimation of the vertical op- 
tical depth Tx(r) = K^Lif). The mean R-band optical depth of 
AUMic is in the range 10" 3 to 10 4 (Kalas 120041 ; using Eq.[3] 
and an area-weighting of r, one obtains an estimate of the ex- 
tinction coefficient (assumed to be constant) 



0.5 M e 



and a final R-band optical depth distribution 



(7) 



(8) 



For a typical inner disk radius of 0.03 AU, t r is smaller than 
unity even in the disk center. A thin disk with such low opti- 
cal depths produces relative deviations of the light curve of the 
order of 0.1%. Current ground-based microlensing observations 
require an order of magnitude higher deviation and so could de- 
tect an average optical depth > 5 • 10~ 3 ; according to Augereau & 
Beust (2006), this would be the maximal optical depth in the vis- 
ible for the AU Mic debris disk, transversed in vertical direction, 
but for an edge-on orientation, optical depths of 4 • 10~ 2 could be 
reached. Consequently edge-on configurations are much more 
likely to be detected from microlensing, too. 

3.2. Geometrically thick disks 

The attenuating effect of a circumstellar disk depends especially 
on the path length of the light ray through the disk and thus on its 
shape. For an edge-on configuration and a 0.5 M star of 0.6 R e , 
the deflection causes a vertical shift of 10~ 3 AU at 150 AU in the 
vertical direction and can be neglected. In the following, a ge- 
ometrically thick disk with a gaussian vertical dust distribution 
and a scale height H(r) = 0.47 (5^7) + 0.34 AU is used (for AU 
Mic cf. Metchev et al. l2005l and Krist et al. 2005). The density p 
of the disk is given by 



p(r,z)= m — g-^W 
ff(r)V2lr 



(9) 



The normalization constant is determined by integrating p paral- 
lel to the disk normal. 



t r = 5.07 ■ 10" 



p(r, z)ds 



(10) 



In PSPL models the light path is usually parameterized using two 
straight lines according to the thin lens approximation: To first 
order, the true light path is a hyperbola: according to Gonner 
( 1 19961 p. 206), the maximal deviation of the distance of closest 
approach R from the hyperbola is 



AR 



2GMR 
Rc 2 + 2GM 



where M is the mass of the deflecting star. The maximal devi- 
ation from the thin lens approximation for grazing incidence is 
10~ 5 AU at 0.6 R e and a deflector mass of 0.5 M Q . 

The sensitivity of different models for r ~ 5 ■ 10~ 4 is shown 
in Fig. [4] The source track for each pair of (z, (f>) was limited to 
a range from -1 to 1 0g and the impact parameter uq = 0.1 6*g 
has been kept fixed; it does not change the maximal variation 
of the transmissivity, only the maximal distance from the cen- 
ter. While the two dimensional model (projected surface den- 
sity) with S oc r -0-8 shows a maximal deviation below 1 %, the 
more realistic three dimensional model with S oc r~° 8 shows an 
increased deviation at high inclinations, because the transversed 
matter is taken into account. The smaller the outer radius of the 
disk is, the poorer the detectability. From the calculated sensi- 
tivity pattern the fraction of the parameter space exceeding 1%, 
i.e. lOmmag, can be calculated: 17% for ro = 50 AU, 34% for 
r D = 100 AU and 50% for r D = 150 AU. These results ar e valid 
for a system without an inner clearing but Wyatt et al. (2007) 
have concluded from observations that 2 + 2% show hot dust 
close to the star. Assuming an inner radius of 10 AU changes the 
detectability to 10% for r D = 50 AU , to 18% for r D = 100A£/ 
and to 23% for r^ = 150 AU. Changing the scale height redis- 
tributes the sensitivity to smaller inclinations. 

In Fig. [5] five different sensitivity configurations are calcu- 
lated for ro = 5QAU and an inner gap with 10AU radius. 
Simulations for 100 AU and 150 AU show similar results: all 
of them show a systematic decrease in the residuals close to 
the maximal magnification, especially at high inclination. The 
source and image tracks plotted over the map of transmissivity 
explain this behavior for a surface density modeled as an inverse 
power law. The image with positive parity reaches its maximal 
separation from the lens when the source-lens distance is min- 
imal and therefore the optical depth oc Z is also minimal. This 
indicates one possibility of detecting circumstellar disks - fitting 
a PSPL model and analyze the structure of the residuals within a 
few Einstein radii. The simulated light curves show deviations 
between 0.1% and 1% (i.e. 1 and lOmmag without blending 
flux), a signal that is detectable by a dedicated space mission 
like the Microlensing Planet Finder or EUCLID. 

4. Relative extinction curves 

It is possible to tell a priori if a certain microlens is surrounded 
by attenuating matter if a wavelength-dependent extinction sig- 
nal can be seen. When fitting a PSPL model, the residuals can be 
used to show where relative attenuation in the lens plane occurs 
between images taken at different wavelengths. A combination 
of photometric and astrometric follow-up could then infer the 
extinction structure of the central zone by measuring the extinc- 
tion coefficient for each image: the centroid of the lensed source 
star will move slightly to the less dimmed image. The centroid 
C with extinction is located at 



9 C = 



k+Q+n + + 



(12) 



(11) 



k + fi + + k-fi- 

where k ± denote the measured transmission coefficients, fi ± the 
unextincted magnifications, 0± the image positions and 6 C the 
trajectories on the sky. The total magnification is described as 

= k+/i + + k-H- (13) 

One of the transmission coefficients can be calculated from Eqs. 
[T2]and[T3] 

H 0c~0 + 



fi- - 6» + 



(14) 
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Fig. 4. Contour lines of maximal relative deviation from PSPL and a grid of attenuated light curves with mo = 0.1 #£ are plotted for 
a surface density S oc r -0 ' 8 , different outer disk radii ro - 50, 100, 150 AU, an inner disk radius of 10 AU, and a constant mean 
optical depth of f = 5 ■ 10 4 



Ideally brightness and position are measured simultaneously and 
thus the measured magnification /j and position 9 C can be jointly 
used in Eq. [14] after correcting for blended flux and baseline. 
Finally k+ is calculated from Eq.[13]by using the image positions 
± from the fitted PSPL model. 

5. Conclusions 

We have shown that gravitational microlensing can - in prin- 
cipal - be used as a tool for detecting circumstellar disks be- 
yond 1 kpc when the photometric residuals relative to a PSPL 
model are searched for deviations between -3 and 3 #g around 
the maximal magnification. We estimate that 4% of all F, G 
& K-microlenses contain debris disks. However, if the optical 
depth t > 5 ■ 10~ 4 (age ~ 1 Gyr), a tenth of these objects should 
have detectable disks. Assuming a power law model for the sur- 
face density X oc r~° 8 and depending on the inner and outer 
radii distribution of debris disks between 10 and 50% of all re- 
maining light curves could show deviations larger than 1%. For 
~ 4000 existing light curves from the OGLE collaboration, we 
expect that of the order of 10 of the F, G & K-lenses show a 
detectable signature induced by a debris disk. If we take into 
account that anomalous microlensing events, i.e. events that are 
not well-described by a PSPL model, can also display systematic 
variations and blending occurs, it is unlikely, that debris disks 
can be discovered from the existing datasets. Nevertheless, fu- 
ture space missions will almost certainly provide the required 
accuracy for detecting debris disks and probing the attenuating 
matter. 
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Fig. 5. Differential light curves of a geometrically thick disk with a central hole are shown for different inclinations i and orientations 
4> of the source track. The disk has an inner radius of 10 AU and an outer radius of 50 AU, X cc r~° 8 ; H = 2.5 AU, and the average 
optical depth r ~ 5 ■ 10~ 4 



